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INTRODUCTION 


Linear  base  sound  ranging,  using  the  GR-8  or  AN/TNS-10  Sound  Ranging  Set, 
has  been  shown  to  be  an  effective  means  for  locating  enemy  firing  batteries 
which  are  beyond  the  range  of  direct  visual  observation.  Development  of 
the  USRAN3  solution  technique  by  Swingle,  Bellucci,  and  Crenshaw  showed 
that  the  present  system  was  not  performing  to  the  limits  of  its  accuracy 
potential.  Their  evaluations  showed  that  an  improvement  in  accuracy  from 
2.2  percent  of  range*  to  1.8  percent  could  be  achieved  by  a simple  change 
in  the  method  by  which  the  final  target  location  (fix)  was  computed. 1 ♦r 
They  noted  that  further  improvements  were  believed  to  be  achievable. 

In  response  to  a request  from  US  Army  Materiel  Command  Headquarters  in 
June  1974,  Swingle  developed  an  error  analysis  for  the  overall  system 
which  was  used  to  estimate  the  likely  improvement  which  could  be  achieved 
for  a typical  target  located  centrally  to  the  microphone  array  and  10,000 
meters  ahead  thereof.  This  restriction  of  target  location  was  adopted 
for  reasons  of  both  simplicity  and  nonavailability  of  quantitative  error 
sensitivity  data  for  other  locations. 

The  present  report  covers  the  first  steps  in  expanding  that  analysis  to 
include  the  entire  field  of  potential  target  locations,  including  exten- 
sion of  the  analysis  to  points  beyond  the  usual  range  of  sound  ranging 
utility,  i.e  , to  ranges  exceeding  twice  the  base  length  and  to  flanking 
angles  beyond  15  degrees  from  the  normal  of  the  outlying  subbases. 

This  analysis  will  develop  the  ultimate  accuracy  which  can  be  expected 
from  a variety  of  sound  ranging  array  configurations  having  tactical 
practicality,  including,  of  course,  the  six-microphone  "linear  base" 
system  with  which  the  US  Army  is  now  equipped.  Optimized  meteorological 
correction  methods,  which  take  account  of  the  natural  variability  of  the 
atmosphere,  will  be  considered  as  parts  of  the  overall  sound  ranging  system 
input  error  source. 

It  will  be  assumed  that  the  reader  is  generally  acquainted  with  the 
principles  of  sound  ranging  on  artillery  weapons  based  on  the  use  of  the 
AN/TNS-10  or  similar  system.  All  six  microphones  will  be  assumed  to  be 
installed  along  a straight  base  at  uniform  spacings  of  four  sound  seconds 
(about  1350  meters).  While  solutions  for  shot  time  and  shot  coordinates 
are  unique,  within  the  usual  assumptions  of  an  atmosphere  having  uniform 


♦Distance  target  is  forward  of  line  of  microphones 

'Donald  M.  Swingle,  Craig  M.  Crenshaw,  and  Raymond  Bellucci,  1972, 
"Improved  Sound  Ranging  Location  of  Enemy  Artillery,"  Army  Science 
Conference,  US  Military  Academy,  West  Point,  NJ 

2Donald  M.  Swingle  and  Raymond  Bellucci,  1973,  "Improved  Sound  Ranging 
Location  of  Enemy  Artillery,"  R&D  Technical  Report  ECOM-5486,  US  Army 
Electronics  Command,  Fort  Monmouth,  NJ 
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wind  and  uniform  temperature  (sound  speed),  when  arrival  times  from  only 
three  microphones  are  used  to  determine  target  location,  the  system  is 
over-determined  when  more  than  three  arrival  times  are  used.  Under  these 
conditions  a number  of  solution  methods  can  be  used,  yielding  generally 
different  target  locations. 

Following  the  development  and  evaluation  of  the  USRAN3  solution,  Crenshaw 
and  Swingle  engaged  in  a continuing  verbal  and  correspondence  exchange,3*4 
looking  forward  to  a more  general  and  even  more  accurate  solution  and  mete- 
orological correction  method.  This  led  to  the  evolution  of  the  Weighted 
Range  Equation  Solution  (WRAS)  which  is  outlined  in  the  appendix.  In  brief, 
this  method  determines  the  shot  time  and  coordinates  which  minimize  the 
sum  of  the  squares  of  the  differences  between  the  square  of  the  geometric 
range  from  solution  point  to  each  microphone  and  the  square  of  the  "prop- 
agation range"  computed  from  the  speed  of  sound  along  that  direction  and 
the  difference  between  computed  shot  time  and  measured  arrival  time. 

N 

W.  x (Pi  “ pf)  = m3ni’mum  (1) 

i=l  ' ' 


where  is  a weighting  factor  which  may  be  varied  by  any  desired  rule 

across  the  microphone  array.  A simple  form  for  the  Wi  which  was  suggested 

by  the  known  behavior  of  atmospheric  variability  as  studied  by  Arnold  and 
Bellucci5  and  later  by  Lowenthal  and  Bellucci6  is 


(2) 


Personal  communications  between  Dr.  Craig  M.  Crenshaw  and  Dr.  Donald 
M.  Swingle,  June  1972-December  1973 

4Craig  M.  Crenshaw,  1972,  "Sound  Ranging  Calculations,"  AMC  Chief 
Scientist  Technical  Note  Number  1,  US  Army  Materiel  Command,  Alexandria,  VA 

5A.  Arnold  and  R.  Bellucci,  1957,  "Variability  of  Ballistic  Meteoro- 
logical I'urameters,"  Technical  Memorandum  M-1913,  US  Army  Signal  Engineer- 
ing Laboratories,  Fort  Monmouth,  NJ 

6Marvin  J.  Lowenthal  and  Raymond  Bellucci,  1970,  "Variability  of 
Ballistic  Winds,  R&D  Technical  Report  ECOM-3259,  US  Army  Electronics 
Command,  Fort  Monmouth,  NJ 
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where  n Is  any  real  number.  Experimentation  was  conducted  with  n ranging 
from  0 to  -10.  When  applied  to  actual  sound  ranging  data,  n = -5  generally 
gave  the  most  accurate  target  locations.  In  the  present  study  this  weight- 
ing was  used,  but  the  conclusions  are  equally  valid  when  solutions  are 
computed  using  n = 0, 


Earlier  studies  of  sound  ranging  system  error  sensitivity  had  been  con- 
ducted by  Fox,7  Lee,8  and  Bellucci ,9»10 

All  the  above  utilized  solutions  depended  upon  use  of  differences  in 
arrival  time  at  adjacent  microphones.  Fox  used  a least  squares  minimiza- 
tion of  the  differences  between  "true"  and  perturbed  arrival  time  differ- 
ences, while  Lee  minimized  the  squared  differences  between  "true"  and 
perturbed  arrival  times.  Bellucci  used  the  standard  artillery  solution. 
Bellucci's  published  work  was  confined  to  a set  of  only  12  points  in  the 
target  area.  Lee  expressed  errors  as  radial  errors  (i.e.,  miss-distance), 
thus  precluding  consideration  of  the  vector  nature  of  the  error  patterns. 

Fox  considered  radial  and  tangential  errors  relative  to  the  center  of  the 
microphone  array.  Lee  presented  total  radial  error  contours  but  provided 
no  data  on  vector  error.  None  of  the  researchers  provided  a dense, 
ordered  array  of  errors  at  target  points  throughout  the  battle  area  which 
might  be  utilized  for  error  analysis  of  both  current  and  future  tactical 
sound  ranging  systems. 

Results  obtained  herein  are  in  generally  good  agreement  with  these  studies. 
They  differ  in  that  they  present  error  transfer  properties  of  the  WRAS 
solution  in  terms  of  along-array  and  along-normal  components.  They  also 
extend  the  range  of  the  computations  to  at  least  five  sound  base  (25  subbase) 
lengths  ahead  of  the  array  and  to  9 to  18  subbase  lengths  to  either  side  of 
the  centerline  of  the  array.  Additional  computations,  not  covered  in  de- 
tail here,  extend  the  computations  out  as  far  as  15  sound  base  (75  subbase) 
lengths  forward  of  the  array.  The  properties  of  linearity  and  superpos- 
ability  discussed  below  were  also  found  to  extend  to  such  ranges,  although 


7H.  L.  Fox,  1968,  "Meteorological  Techniques  for  Sound  Ranging;  Theory 
of  Errors,"  Technical  Report  EC0M-0233-2,  US  Army  Electronics  Command, 

Fort  Monmouth,  NJ 

8Robert  P.  Lee,  1972,  "Artillery  Sound  Ranging  Computer  Simulations," 

R&D  Technical  Report  ECOM-5441  , US  Army  Electronics  Command,  White  Sands 
Missile  Range,  NM 

9Raymond  Bellucci,  1966,  "Studies  of  Meteorological  Techniques  for 
Sound  Ranging:  Report  II,  Error  Analysis,"  Technical  Report  EC0M-2703, 

US  Army  Electronics  Command,  Fort  Monmouth,  NJ 

10Personal  communications  between  Raymond  Bellucci  and  Dr.  Donald  M.  Swingle, 
June  1966- June  1972 
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errors  were  generally  much  larger.  These  results  will  be  utilized  in  suc- 
ceeding reports  dealing  with  errors  in  microphone  position  and  atmospheric 
input  as  well  as  terrain  effects  and  sound  ranging  system  improvement. 


DATA  PRESENTATION 

i 

The  data  of  this  report  are  uniformly  presented  as  fields  of  numerical 
values,  to  provide,  in  a reasonably  compact  form,  the  detailed  informa- 
tion required  to  assess  the  effects  of  various  error  sources  on  target 
location.  Those  readers  whose  taste  leads  to  the  representation  of  data 
in  contour  maps  are  invited  to  add  their  own  analyses  of  the  numerical 
fields  presented.  Both  horizontal  components  of  target  location  error 
are  presented  at  each  grid  point  to  facilitate  ready  analysis  of  the 
vector  character  of  the  target  location  error  associated  with  input  data 
errors . 


| DATA  INTERPRETATION 

Figure  1 provides  a schematic  guide  to  interpretation  of  the  following 
tables.  Tables  1 through  23  provide  the  field  of  target  location  errors 
due  to  specified  errors  in  arrival  time  at  given  microphones  for  each  of 
a regular  array  of  assumed  "true"  target  locations.  The  targets  are 
assumed  to  be  at  the  points  indicated  by  "+"  (except  along  the  left-hand 
edge  of  the  figure). 


These  points  are  uniformly  spaced  by  one  subbase  length  forward  of  the 
center  of  each  subbase  and  laterally  from  the  center  of  the  array.  To 
the  left  of  each  target  point  appear  two  figures.  The  upper  number  is 
the  error  in  meters  which  would  arise  in  the  solution  by  WRAS  for  the 
target  along  the  X-  or  along  baseline  direction,  while  the  lower  number 
is  the  error  in  the  Y-  or  normal-to-baseline  direction. 

Along  the  top  of  each  figure  are  listed  the  subbase  length  (BL)  in  sound 
seconds  (4  sound  seconds  is  approximately  1350  meters)  and  the  index  number 
of  the  perturbed  microphone  (IP,  counting  from  left  to  right  as  one  faces 
the  target  from  the  baseline).  In  this  report  IP  will  range  from  1 to  6. 

If  two  digits  are  given  under  IP,  these  digits  indicate  the  two  micro- 
phones whose  arrival  times  are  simultaneously  perturbed.  ID  is  used  to 
indicate  the  perturbed  parameter  and  will  always  be  8 in  this  report, 
indicating  that  the  perturbation  considered  is  in  time  of  arrival  at  the 
microphones  indicated  by  IP.  The  number  of  subbase  lengths  to  the  right 
of  the  array  center  which  characterizes  the  top  row  of  targets  (+)  is 
given  by  M.  As  the  computations  progress  down  the  page,  M is  incremented 
for  each  row  and  appears  at  the  left-hand  edge  of  the  paper  abreast  of  the 
Y-component  of  error.  Just  above  these  changing  values  of  M appears  the 
value  of  IP  given  in  the  heading  line.  This  value  appears  just  to  the 
left  of  the  first  + in  each  row  which  marks  the  point  which  is  M subbase 
lengths  from  the  array  center.  In  all  of  the  tables  N = 25,  indicating 
that  the  right-most  + or  target  position  in  each  row  is  25  subbase  lengths 
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forward  of  the  microphone  array.  DXY  indicates  the  spacing  in  both  X and 
Y directions  of  the  grid  points  in  units  one  subbase  long.  In  all  tables 
herein  DXY  is  1.000.  UP  gives  the  assumed  unit  perturbation  which  is 
applied  to  the  arrival  time  of  the  identified  microphone(s)  expressed  in 
milliseconds.  IW  is  the  weighting  exponent.  In  terms  of  equation  (2), 
n = -IW,  IW  will  always  be  5 herein.  ZZ  gives  the  target  height  above 
the  plane  containing  the  microphones  and  is  always  0.0  in  this  report. 

VI,  V 2,  and  V3  are  the  X-,  Y-,  and  Z-  components  of  the  mean  wind  assumed 
in  the  computation,  and  will  remain  0.00  in  this  report.  Finally,  TE  is 
the  assumed  "effective  temperature"  for  sound  ranging,  assumed  to  be  10°C 
herein. 

Where  ****  appears  on  the  tables,  the  error  exceeds  ±999.5  meters,  while 

"q+  appearing  in  the  upper  and  lower  left  corners  indicates  that  computa- 
tions were  not  made  for  the  point  because  it  occurs  at  a large  flanking 
angle  ( | M J > N + 3) . 


EFFECT  OF  UNIT  INPUT  ERRORS  ON  SOLUTIONS 

In  computing  the  data  for  tables  1 through  15,  the  author  assumed  that 
the  target  was  at  the  plotted  grid  point,  and  the  arrival  time  of  sound 
generated  at  shot  time  (=  0,000  seconds)  was  computed  for  each  microphone. 
Then  the  arrival  time  at  the  specified  microphone  (IP)  was  incremented  by 
UP  milliseconds,  and  the  target  location  was  computed  by  using  the  WRAS 
solution.  At  the  field  point  the  difference  between  the  computed  location 
and  the  assumed  or  "true"  location  is  printed  (computed  - true).  In  tables 
1 through  6 of  "error  maps,"  the  unit  perturbation  was  assumed  to  be  1 mil- 
lisecond, All  errors  are  rounded  to  the  nearest  meter.  This  value  was 
chosen  as  small  enough  to  produce  few  errors  exceeding  1000  meters,  while 
also  producing  few  errors  less  than  a few  meters  in  at  least  one  of  the 
component  directions.  The  reader  is  invited  to  examine  the  symmetry  of 
error  effects  between  microphones  1 and  6 (tables  1 and  6),  2 and  5 
(tables  2 and  5)  and  3 and  4 (tables  3 and  4). 

The  tables  show  that  errors  in  arrival  time  at  microphones  2 and  5 have 
relatively  small  effects  at  most  target  points  on  either  component  of 
target  location  error  and  that  for  all  microphones  the  Y-component  error 
is  generally  several  times  the  X-component  error.  Further,  both  micro- 
phones 1 and  6 yield  large  negative  Y-component  errors  for  a small  posi- 
tive change  in  arrival  time,  while  both  microphones  3 and  4 yield  large 
positive  errors  in  Y-component  for  a small  positive  change  in  arrival 
time. 

In  the  above  computations  and  throughout  this  report,  it  does  not  matter 
whether  the  "error"  in  arrival  time  be  the  result  of  an  error  in  reading 
the  acoustic  record  or  the  effect  of  propagation  phenomena  in  causing 
the  signature's  apparent  arrival  time  to  be  different  from  that  which 
would  be  expected  after  correcting  for  the  assumed  meteorological  condi- 
tions. Either  cause  will  result  in  the  same  target  location  error  if  it 


causes  the  same  error  in  arrival  time.  For  example,  at  the  point  M = -9, 
N ■ 25,  the  upper  right-hand  + in  table  1,  one  finds  that  an  arrival  time 
error  of  1 millisecond  (UP  = 1„00)  causes  the  X-component  of  target  posi- 
tion to  be  computed  to  be  18  meters  larger  than  true  and  the  Y-component 
to  be  computed  to  be  49  meters  smaller  than  true  (18,  -49).  For  the  same 
point,  tables  2 through  6 show  that  the  same  unit  error,  if  it  occurs  at 
microphones  2 through  6,  produces  target  location  errors  of  (-5,  17), 

(-14,  40),  (-11,  31),  (-1,  1),  and  (13,  -41)  meters,  respectively. 


LINEARITY  OF  ERROR  TRANSFER  FUNCTIONS 

The  analysis  of  the  joint  effects  of  several  error  sources  on  the  overall 
sound  ranging  target  location  error  would  be  much  simplified  if  the  output 
error  could  be  assumed  to  be  linearly  dependent  upon  the  input  error,  or  at 
least  approximately  so.  Tables  2 through  6 show  the  error  fields  for  input 
* errors  of  1 millisecond  occurring  in  microphones  1 through  6,  while  tables 

7 through  12  show  the  error  fields  for  input  errors  of  3 milliseconds. 

I 

In  view  of  the  symmetry  noted  for  errors  in  microphones  1 and  6,  2 and  5, 
and  3 and  4,  it  becomes  apparent  that  the  entire  story  can  be  told  by  just 
j the  error  fields  due  to  microphones  1,  2,  and  3.  Thus  the  error  fields 

due  to  input  errors  of  10  and  30  milliseconds,  shown  in  tables  13  through 
15  and  16  through  18,  respectively,  are  shown  only  for  microphones  1,  2, 
and  3. 

Through  most  of  the  error  field,  output  errors  in  each  horizontal  component 
of  target  location  are  approximately  proportional  to  input  errors.  Over  a 
range  of  30:1  in  input  error,  proportionality  is  quite  good.  The  range 
selected  for  illustration  includes  the  range  of  errors  to  be  expected 
because  of  normal  atmospheric  propagation  effects  and  from  errors  in 
reading  arrival  time  records. 

In  assessing  the  combined  effects  of  several  errors  affecting  the  measured 
arrival  time  at  any  microphone,  little  error  will  occur  if  the  target 
location  error  due  to  each  contributing  error  is  separately  computed  and 
the  results  summed  instead  of  first  summing  the  several  contributing 
errors  in  arrival  time  and  then  computing  the  target  location  error  which 
they  will  cause. 


SUPERPOSABILITY  OF  ERROR  TRANSFER  FUNCTIONS 

The  preceding  section  demonstrated  that  several  simultaneous  arrival  time 
errors  occurring  at  a single  microphone  could  be  separately  transformed 
into  target  location  errors  and  then  summed  to  give  the  aggregate  target 
location  error.  To  investigate  the  validity  of  computing  errors  due  to 
simultaneous  arrival  time  errors  in  several  microphones  as  the  sum  of  the 
effects  of  errors  in  each  microphone,  the  author  assumed  errors  occurring 
in  one  microphone  and  an  additional  error  occurring  in  each  of  the  remain- 
ing microphones  in  turn. 
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In  tables  19  through  23,  simultaneous  errors  of  3 milliseconds  are  assumed 
to  occur  in  microphone  2 and  in  microphones  1,  3,  4,  5,  and  6 in  turn. 

The  reader  will  be  able  to  readily  verify  that  the  effect  of  the  combined 
simultaneously  occurring  error  is  very  closely  approximated  by  the  sum  of 
the  errors  given  table  8 and  tables  7,  9,  10,  11,  and  12,  respectively. 


For  example,  tables  7 through  12  show  that  the  errors  in  target  location 
at  the  point  (M,  N)  = (-9,  25)  due  to  a 3-millisecond  error  in  arrival 
time  at  microphones  1 through  6 are  (56,  -147),  (-15,  50),  (-42,  120), 

(-34,  93),  (-3,  2),  and  (40,  -124)  meters,  respectively.  For  the  same 
target  point,  table  19  shows  that  the  target  location  error  due  to  simulta- 
neously occurring  3-millisecond  errors  in  microphones  1 and  2 is  (40,  -95), 
almost  exactly  the  sum  of  the  effects  of  this  error  acting  separately  at 
microphones  1 and  2.  The  corresponding  effect  of  the  combined  error  in 
microphones  2 and  3,  2 and  4,  2 and  5,  and  2 and  6 are  (-58,  172), 

(-50,  145),  (-19,  53)  and  (15,  -74)  meters,  respectively,  as  shown  in 
tables  20  through  23.  In  each  case  the  effect  of  the  simultaneously 
occurring  error  is  approximately  equal  to  the  sum  of  the  effects  caused 
by  each  error  taken  as  acting  alone.  Similar  results  would  be  found  for 
any  selected  group  of  microphones  and  for  any  target  point. 


We  have  thus  shown  that  the  errors  in  target  location  resulting  from  the 
simultaneous  presence  of  errors  in  several  microphones  can  be  closely 
approximated  by  the  sum  of  errors  computed  as  if  each  error  occurred 
separately  in  each  of  the  microphones. 


CONCLUSIONS 


The  above  analysis  has  shown  that  sound  ranging  solutions  obtained  by  the 
WRAS  solution  technique  applied  to  a regular,  linear,  six-microphone  array 
have  the  following  properties: 


1.  Errors  occurring  in  arrival  time,  relative  to  the  assumptions 
inherent  in  the  solution  model,  in  microphones  1,  3,  4,  and  6 have  the 
greatest  effect  on  the  solution;  while  at  most  field  points,  microphones 
2 and  5 have  relatively  little  effect  on  the  solution. 


2.  At  most  field  points,  the  Y-component  of  target  location  error 
due  to  an  error  in  arrival  time  at  any  microphone  is  generally  several 
times  as  large  as  the  X-component  of  error. 


3.  Positive  errors  in  arrival  time  at  microphones  1 and  6 lead 
generally  to  negative  Y-component  errors,  while  positive  errors  in  arrival 
time  at  microphones  3 and  4 result  in  positive  Y-component  errors. 


4.  Positive  errors  in  arrival  time  at  microphones  1 and  6 lead  to 
positive  X-component  errors  in  target  locations  for  left  flank  targets 
and  negative  X-component  errors  in  target  locations  for  right  flank  tar- 
gets, while  positive  errors  in  arrival  time  at  microphones  3 and  4 lead 
to  positive  X-component  errors  in  target  location  for  right  flank  targets 
and  negative  X-component  errors  in  target  location  for  left  flank  targets. 
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5,  For  each  horizontal  component  of  target  location,  the  error  caused 
by  an  error  in  arrival  time  at  any  microphone  is  approximately  proportional 
to  the  arrival  time  error;  i.e.,  the  time-to-space  error  transfer  function 
ts  approximately  linear.  Thus  the  error  in  target  location  due  to  the  com- 
bined effects  of  several  errors  in  arrival  time  at  any  one  microphone  is 
closely  approximated  by  the  sum  of  errors  in  target  location  which  each 

of  the  contributing  errors  would  have  caused  if  acting  alone. 

6.  The  error  in  target  location  in  each  component  of  position  result- 
ing from  simultaneously  occurring  errors  at  two  or  more  microphones  is 
closely  approximated  by  the  sum  of  the  errors  which  each  arrival  time 
error  would  produce  if  acting  alone;  i.e.,  the  error  contributions  are 
superposable  in  geometric  space. 


TABLE  1.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 1-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  1 


TABLE  3.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 1- MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  3 


TABLE  4.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 1-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  4 


TABLE  5.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 1- MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  5 


TABLE  6.  X-  ANO  Y-COMPONENT  ERRORS  DUE  TO  A 1- MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  6 
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TABLE  7.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 3-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  1 


TABLE  8.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 3-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  2 


TABLE  9.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 3-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  3 
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TABLE  10.  X-  AND  Y-C0MP0NENT  ERRORS  DUE  TO  A 3-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  4 
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TABLE  11.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 3-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  “TCROPHONE  5 


TABLE  12.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 3-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  6 


TABLE  13.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 10  MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  1 


TABLE  14.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 10-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  2 
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TABLE  15.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 1 O-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  3 


TABLE  16.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 30-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  1 


TABLE  17.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 30-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  2 


TABLE  18.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  A 30-MILLISECOND  ERROR  IN  ARRIVAL  TIME  AT  MICROPHONE  3 


TABLE  19.  X-  AND  Y-COUPONENT  ERRORS  DUE  TO  3-MILLISECOND  ERRORS  IN  ARRIVAL  TIME  AT  MICROPHONES  1 AND  2 


TABLE  20.  X-  AND  Y-C0MP0NENT  ERRORS  DUE  TO  3-MILLISECOND  ERRORS  IN  ARRIVAL  TIME  AT  MICROPHONES  2 AND  3 


TABLE  21.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  3-MILLISECOND  ERRORS  IN  ARRIVAL  TIME  AT  MICROPHONES  2 AND  4 


TABLE  22.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  3-MILLISECOND  ERRORS  IN  ARRIVAL  TIME  AT  MICROPHONES  2 AND  5 


TABLE  23.  X-  AND  Y-COMPONENT  ERRORS  DUE  TO  3-MILLISECOND  ERRORS  IN  ARRIVAL  TIME  AT  MICROPHONES  2 AND  6 
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APPENDIX 

THE  WEIGHTED  RANGE  EQUATION  SOLUTION 


Given  microphone  positions  (x,  y,  z).,  arrival  time  t^ , target  position 
(X,  Y,  Z),  shot  time  T,  effective  mean  air  velocity  componets  (u,  v,  w)i 
and  sound  speed  c. , we  can  write  the  equation  connecting  these  as  follows: 

{(X  - x,)  ♦ u.(ti  - T)J2  + [(Y  - y,)  * v,(t,  - T)] 

* [(z  - zi)  * "i(‘i  - t)]2  - cK(i  • T)2  = et  • (fll) 


With  real  data,  there  may  be  errors  in  our  determination  of  any  of  the  sub- 
scripted variables  in  the  left-hand  member,  so  that  the  e.  are  generally 

nonzeroc  Expansion  and  rearrangement  yield 


R?  - P?  = e.  , 
i i 1 


where,  denoting  the  effective  speed  of  sound  travel  toward  each  microphone 
by  sit 

• (*  - *02  + (Y  - *i)2  * (z  - 2t)2  • <A3> 

Pi  = si(ti  - T)2  » (A4) 

s!  ' (ct  - u)  - vi  - wi  ) ■ 2/(li  - T) 

x [u/x  - x.)  + v. (y  - y,)  + wi(z  - i,)  • <A5) 
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Thus  defined,  the  weighted  range  equation  solution  seeks  those  values  of 
the  unknown  target  coordinates  and  shot  time  (X,  Y,  Z,  T)  for  which 


N 

y W.j e?  = minimum  » (A6) 

i=l 

For  the  purposes  of  this  study,  y.. , , u^ , , and  Z are  taken  to  be 

zero,  ail  c^  are  assumed  equal  to  a single  value  (the  speed  of  sound  at  an 

effective  temperature  of  10°C),  and  x^+^-xi  is  assumed  to  be  a constant, 

1350  meters,  Based  on  prior  experimentation  with  real  sound  ranging  data 
to  determine  the  best  exponent  n,  W^.  is  defined  as  follows: 

/ N 

wi  = R?  / Z R-  » (A7) 

/ i=1 

n = -5  , (A8) 


The  solution  to  (A6)  is  obtained  by  an  iterative  computation,  using  con- 
ventional least  squares  methods  and  treating  the  W.  as  constants  during 
each  approximation,  1 
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